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Trefftz discontinuous Galerkin methods 
on unstructured meshes for the wave equation 
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Abstract — We describe and analyse a space-time Trefftz discontinuous Galerkin method for the wave equation. The method is 
defined for unstructured meshes whose internal faces need not be aligned to the space-time axes. We show that the scheme is 
well-posed and dissipative, and we prove a priori error bounds for general Trefftz discrete spaces. A concrete discretisation can 
be obtained using piecewise polynomials that satisfy the wave equation elementwise. 
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1 Introduction 

Finite element-type schemes whose discrete functions are 
piecewise solutions of the discretised PDE are often named Tr¬ 
efftz methods. When applied to time-dependent problems, Tre¬ 
fftz methods require the use of a space-time mesh, as opposed 
to the combination of a discretisation in space and a time¬ 
stepping. In the context of linear wave propagation, Trefftz 
methods have been widely studied in time-harmonic regime 
(e.g. |4|), while only few recent works have been devoted to 
the time-domain acoustic and electromagnetic wave equations, 
chiefly m 0 0 0 ED- Some possible advantages of Tre¬ 
fftz methods over standard ones are better orders of conver¬ 
gence; flexibility in the choice of basis functions; low disper¬ 
sion; incorporation of wave propagation directions in the dis¬ 
crete space; adaptivity and local space-time mesh refinement. 

Here we introduce a Trefftz discontinuous Galerkin (DG) 
method for the wave equation written as a first-order system. 
We consider unstructured meshes whose faces are not aligned 
to the space-time axes. The metod proposed is an extension of 
those in mm to higher dimensions and more general meshes, 
and is closely related to those in 0 (Maxwell’s equations) and 
in 0 (non-Trefftz, see Remark U below). After introducing 
the initial boundary value problem in j|2] some notation in $3] 
and the Trefftz DG formulation in $4] in 45. 2 1 we show that the 
Trefftz-DG formulation is well-posed for any Trefftz discrete 
space and we prove a priori error bounds in a DG norm. In 
45.31 - jl5.5l we study the stability of an auxiliary problem in or¬ 
der to give a sufficient condition for an a priori error bound on 
the space-time L 2 norm of the Trefftz-DG error, and show that 


for a class of meshes this condition is verified. These meshes 
can be obtained with the “tent-pitching algorithm” of mmm, 
and allow to treat the Trefftz-DG method as a semi-explicit 
scheme and to compute the solution element by element by 
solving a sequence of small local linear systems. We also show 
that the proposed method is dissipative. The analysis is carried 
out closely following that of 00; the duality argument used 
to control the L 2 norm of the error follows an idea of J8). 


2 The initial boundary value problem 


We consider an initial boundary value problem (IBVP) posed 
on a space-time domain Q = fl x (0,T), where G C R” is 
a bounded Lipschitz polytope with outward unit normal iiq, 
and T > 0. The boundary of G is divided in three parts, 
denoted T^, and l’ K , corresponding to Dirichlet, Neumann 
and Robin boundary conditions; one or two of them may be 
empty. The wave equation IB VP reads as 


(1) 


_ da 

Vv + — = 0 
dt 


V • a 


1 dv 

Jdi 


= 0 


v(-,0)=v o , <x(-,0) 


v =9 d , 
a ' n n = 9n, 

-v-a ■ n n = g R , 
c 


in Q , 

in Q, 

a o on G, 

on Td x [0,T], 
on Tat x [0,T], 

on T fl x [0, T\. 


Here i’o,a o,gD,9N,9R are the problem data; c > 0 is the 
wave speed, which is assumed to be piecewise constant and 
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independent of f; t? £ L°°(T R x [0,T]) is an impedance pa¬ 
rameter, which is assumed to be uniformly positive. The gra¬ 
dient V and divergence V- operators are meant in the space 
variable x only. The equations in <(T]) may be derived from the 
second-order scalar wave equation — AU + c~ 2 j^pU = 0 set¬ 
ting v = QT- and <t = — TU. 


3 Space-time mesh and notation 

We impose a finite element mesh Th on the space-time domain 
Q. We assume that all its elements are Lipschitz polytopes, so 
that each internal face F = dK\ r\dK 2 , for K\. K 2 £ Th. with 
positive n-dimensional measure, is a subset of a hyperplane: 

Fcn f - {(x,t) 6K" +1 : x-n£ + fn‘ F =C F }, 

where (np,rip) is a unit vector in R” +1 and Cp £ K. We 
assume that all the internal mesh faces belong to one of the 
following two classes: 


on an internal face F = <)K\ n OK 2 , either 

{ c|nfi| < iflp and F is called “space-like” face, or 
iflp = 0 and F is called “time-like” face. 

(See Remark|T|for the extension to more general meshes.) On 
space-like faces, by convention, we choose rip > 0, i.e. the 
normal unit vector (n^., rip) points towards future time. Intu¬ 
itively, space-like faces are hypersurfaces lying under the char¬ 
acteristic cones and on which initial conditions might be im¬ 
posed, while time-like faces are propagations in time of the 
faces of a mesh in space only. We use the following notation: 


Th '■= ^ dl\ (the mesh skeleton), 

J^ pace := the union of the internal space-like faces, 

T l ™ c := the union of the internal time-like faces, 
T° h :=Qx {0}, Tl := Q x {T}, 

T? := r D x [0, T\, .Ff := T N x [0, T), 

T*:= T R x [0, T}. 

(We will consider more specific meshes with T)) me = 0 in 35.51 
below.) We denote the outward-pointing unit normal vector on 
dK by {n!p, rip). For piecewise-continuous scalar (w) and 
vector (r) fields, we define averages ({•]}, space normal jumps 
[■]n and time (full) jumps [-]t in the standard DG notation: on 
F = dK ! n dK 2 , K\, K 2 £ %, 


M 

W 

l*i 

-\-W\ 

2 ’ 

Hn 

:= w 

K 1 

n x 

U K 1 

+ w \K 2 n K 2 i 

Hn 

~T 


' n ifi + t \k 2 ' 11 

M t 

:= w 

Kl 

n± K x 

+ W \ K 2 n K 2 

Hi 

:=t 

*1 

n ifi 

+ T \K 2 n K 2 


\k ± + t \k 2 


n K 2 = ( w - w + )np , 
= ( r~ - r + )n^, 


where, on space-like faces, w~ and w + denote the traces of 
the function w from the adjacent elements at lower and higher 
times, respectively (and similarly for r ± ). We use the nota¬ 
tion |-]n to recall that [w]n is a vector field parallel to nf. and 
[r] N is the jump of the normal component of t only. We re¬ 
call also that in this formulas \n t K \, |n^| < 1, and that one of 
the two might be zero; in particular [w| t = 0 and [r] t = 0 on 
Tl me . The following identities can easily be shown: 

- \W\t = 

F 

(3) r-[r]^i[r 2 ] t = ^-|[r] t | 2 , 

W“[t]n + T ~ • [ui]n - Hn = 4N*Mn, 

Up 

+ {{t}} • [w] N = [wt]n, 

We assume that the mesh Th is chosen so that the wave speed c 
is constant in each element; as c is independent of time, it may 
jump only across faces in T )J me . 

Finally, we define the local and the global Trefftz spaces: 

T (K) :={(«;, r) £ L 2 (K) 1+n , s.t. r {gK £ L 2 (dK) n , 

V • r £ L 2 (K), Vw £ L 2 (K) n , 
v„, + ^ = °.v. r + .Lgg = 0 } VA £ 71, 

-J 

v- 

(Trefftz property) 

T {%) :={(t»,r) £L 2 (Q) 1+n , s.t. 

(w\ K ,T\ K )£T(K)VK£T h }. 

The solution (v , er) of IB VP ([TJ is assumed to belong to T (Th ) ■ 


4 The Trefftz-discontinuous Galerkin method 


To obtain the DG formulation, we multiply the first two equa¬ 
tions of ([]} with test fields r and w and integrate by parts on a 
single mesh element K £ Th- 


<4) “X(”( v r+ ^) +,T (lt +v ”')) <iy 

+ J T + a w) ■ n x K + (cr ■ T + dS = 0. 


We look for a discrete solution (vh p ,crhp) approximating 
(v. rr) in a finite-dimensional (arbitrary, at this stage) Trefftz 
space V(7ji) C T(7 h). We take the test field (w,t) in the 
same space V (Th), thus the volume integral over K in £j} van¬ 
ishes. The traces of Vh p and <T/, r , on the mesh skeleton are 
approximated by the (single-valued) numerical fluxes Vh P and 
a’hp, so that (Q} is rewritten as: 


(5) 




t ■ n 


X 

K 


+ ^2 T crhp • {wrTp + t n^p) dS" — 0. 
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We choose to define the numerical fluxes as: 


Vhp ’ — 

' v hp 

Vhp 

VO 

9d 

v h p+P{(Thp ■ n n~g N ) 

(1 - S)v hp 

. +T (^hp-n^+gn) 


G hp • — 

_ 

a hp 
G hp 

Co 

IM + “Mn 

<j h p+ot(v h p-gD)T 

9n 

k ~\~8cThp 


T-space 

on Fh , 

on J-J, 
on J-o, 
on F™ e , 
onT°i 

onJ-f, 

onJ?. 


The mesh-dependent parameters a £ L 00 (F% me U Fj?), /3 £ 
L°°(F™ e U T^) and 6 £ L°°{T^) may be used to tune the 
method, e.g. to deal with locally-refined meshes. 


The fluxes are consistent, in the sense that they coincide with 
the exact solution (v, cr ) of the IB VP ([U if they are applied to 
( v,cr ) itself, which satisfies the boundary conditions and has 
no jumps across mesh faces. Moreover, the fluxes satisfy the 
boundary conditions, e.g. -Vh P — dh P ■ = gn on F^. The 

numerical fluxes can be understood as upwind fluxes on the 
space-like faces and centred fluxes with jump penalisation on 
the time-like ones. 


Summing the elemental TDG equation (0 over the elements 
K £ Th, with the fluxes defined above, we obtain the Trefftz- 
DG variational formulation displayed in box 0. 

Method © appears as a formulation over the whole space-time 
domain Q, which would lead to an unconvenient global linear 
system coupling all elements. However, if the mesh is suitably 
designed, the solution can be computed by solving a sequence 
of smaller local problem. A first possibility is to partition the 
time interval [0, T] into subintervals and solve sequentially for 
the corresponding time-slabs f 1 x (tj-i,tj), see 0] ! 6); this 
corresponds to an implicit method and allows local mesh re¬ 
finement. A slightly more complicated, but potentially much 
more efficient version is to solve for small patches of elements, 
localised in space and time, in the spirit of the semi-explicit 
“tent-pitching” algorithm of (2] (3j |7j| ■ If no time-like faces are 
present in the mesh, the solution can be found by solving a 
small system for each element, see £ 15. 51 below. If the same 
mesh is used, all these variants are equivalent, in the sense that 
the discrete solutions ( Vh pi crhp) coincide. 


REMARK 1 One could easily extend the formulation weaken¬ 
ing assumption (0 to allow more general time-like faces with 
c|n|,| > rip, i.e. not aligned to the time-axis. Choosing the 
numerical fluxes as above, one obtains a formulation similar to 
0 with additional terms on F™ e containing [ic ] t and [/r]t, 
which do not vanish in this setting. It is then easy to prove 
the coercivity of the new bilinear form in a slight modification 
of the DG norm introduced below. However, the bilinear form 
will contain the term • [r] t dS, featuring the full 


jump of t (as opposed to the normal jump only), which, in di¬ 
mension n > 1, is not easily controlled by the same DG norms. 

Remark 2 Formulation @ can be seen in the framework 
of DG methods for general first-order hyperbolic systems 
developed in (7), which considers standard discontinuous 
piecewise-polynomial spaces. The choice of the numeri¬ 
cal fluxes in the interior elements correspond to the choice 
of a suitable decomposition of the block matrix M = 
,n K c (n K ) \ defined on Qx for all K £ Th- Here Id n 

v n f K nf K Id„ > inn 

is the identity matrix in B, rLXn , T denotes vector transposi¬ 
tion, and n x K is thought as a column vector. The choice we 
have done in this section corresponds to the decomposition: 


on 

M+ = 

M- = 

dK n 

M 

0 

9kn7; pace n{<<0} 
dK n F™ e 

0 

( “ 5 ( n ?c) T \ 

K \*K /3n k®*k ’ 

M 

, —a i(n^) T , 
^ n K —P n K® n K 


The conditions M + + M = M, ker(M + — M ) = ker(M) 
and Ml + Ml = Oon c)K\ fl dK 2 are automatically 

satisfied. Moreover, M + > 0 and M - < 0 hold true if and 
only if a/3 > 1/4. The boundary terms in ((6)) and in E §6] 
coincide (apart from a different sign convention) if our flux pa¬ 
rameters and their boundary coefficients Q and a are chosen 
so that a = cr on Td, /3 = 1/tr on Tjv, 5 = (1 + Q)/2 and 
d/c = cr( 1 + Q)/( 1 - Q) on T fl . 


5 A priori error analysis 
5.1 Definitions and assumptions 


We prove the well-posedness and the stability of the Trefftz- 
DG formulation 0 under the assumption that the flux param¬ 
eters a, j3 and 8 are uniformly positive in their domains of def¬ 
inition and that ||<5|| L00 prR) < 1. We introduce a piecewise- 
constant function 7 defined on J r/ pilcc , measuring how close to 
characteristic cones the space-like mesh faces are: 


(7) 7:= 


cm 


onTc F^ ce , 


7 := 0 onJ/U F£, 


from which, recalling assumption 0, 7 £ (0,1) and 
(8) IMnI < ~ c IHtl, |[t]n| < ~ I[t] 1 1 onJT ce . 

We define two mesh- and flux-dependent norms on T(7jj): 

2 


HGT)||| dg := 

■1 -7\V2 


T — 7\ !/2 , 

' C >] t 


(¥) 


L 2 (JF/ pa “) 


(¥)' 


Fit 


-1 11 2 

2 ll c w \\L 2 (r° u ^) ^ 2 11 ' 


L 2 (J r / pace ) n 
1 11 11 2 

+ ~ M 


o 1/2 Hn 






N 


2 

L 2 (^ me ) 
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+ 


(A) 


1/2 


L 2 (^ pa “)” 


- Mn - ~[c 2 w 2 + |r| 2 ] t ) dx 


+ 

r 1/2 w 

2 

+ 

a~ 1/2 irj 



L 2 (JF““) 



+ 

a~ 1/2 r • n£ 

2 

+ 

f)- x/2 w 






2 

L 2 (J^ m *) n 


(The factors (1 — 7) ±1 ^ 2 may be dropped from all terms in 
the two norms without modifying the analysis in the follow¬ 
ing sections.) We note that these are only seminorms on bro¬ 
ken Sobolev spaces defined on the mesh 7/t, but are norms on 
T(7h): indeed |||(u),r)||| = 0 for (w,t) G T(%) implies 
that (w, t) is solution of the IBVP 0 with zero initial and 
boundary conditions, so (w, t) = ( 0 , 0 ) by the well-posedness 
of the IBVP itself (see 0 Lemma 4.1]). 

5.2 Well-posedness 

We first note that for all Trefftz field (w, r) G T(7 ~h) 


+ \ [ (c 2 w 2 + |t| 2 ) dx + i [ (c 2 w 2 + |r| 2 ) dx 
2 ■% 2 

+ f (HMn + W-Nn 

+ a|NN| 2 +/3 [t]n - K]nJ dS 
+ [ aw 2 dS+[ p(r • n£)(r • nfJdS 




(l-5)l9 2 SC 
- w + T 

C V 


(r-n£) 2 ) 


d S 


1 / (^r( c ~ 2 l w it + II T ltl 2 ) + tHiWn ) dx 

V ‘-‘Tip n F 

1 11 _i 112 1 || || 2 

+ 2 H c + 2 


o 1/2 Hn 


IdK 


wt ■ n x K + 




i 2 (^r) 


/3 1/2 M N 

2 


L 2 (tF“ me ) 


0 V/\ G %, 

+ 

a^w 

2 

+ 

/3 1/2 t ■ 







which follows from integration by parts and the Trefftz prop¬ 
erty. Subtracting these terms from the bilinear form A, using 
the jump identities ( 0 , the inequalities 0 , the definition of 7 
in 0 . and the weighted Cauchy-Schwarz inequality, we show 


+ 


HMD 

> 


<N 

1 

r-H 

2 

(ScAA x 

K c ) w 

+ 

^ 2 (^ B ) 

V d ) ' " n 

3 IIK'^'OIIIdG 

V(w 

,t) GT(Tfc). 
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Using again the Cauchy-Schwarz inequality, the bounds on the 
jumps © and 7 < 1 , we have the following continuity estimate 
for the bilinear form A: for all (v, cr), (w, r) G T(7 h) 


where 


= 0 , 

2max | || VIL-p?) ’ lll=5ll/«>(^«)} 7^ 0- 



«, 0 - 1 IIIdg+II|(w,t)|||dg, 


Note that C c > 2 and that the minimal value C c = 2 is ob¬ 
tained for <5 = 1/2. Also the linear functional l is continuous: 

\V(w,t )| <(2 ||c - 1 t;o||^ 2 (:F o ) + 2 ||o- 0 ||^ 2( ^ ) 


+ 2 


a 1/2 g D 


l 2 (f?) 


P 1/2 g. 


(c/d) 1/2 g R 


2 y/ 2 , 

L 2 (r?) 


N 

L 2 (J~hi 
^,t)|||dg+; 

rD ttN 


{V - Vhp,cr - a hp ) |||dg 

< (1 + Cc ) inf I 
(u>,-r)ev(r h ) 


(v — w.cr — t) 


|dg+j 


with C c as in Moreover, if go = gN = 0 then 
11 K^hpj O'Tip) 11 |dG ^ 


(2 

1 

-Vo 

V 

c 


+ 2 ||(To| 


L 2 (F° h ) 


L 2 (K) 


, 1/2 


t ? 1 / 2 




i 2 (^ K ) 


1/2 


( 10 ) 


S7z + dC/dt = <f> 

in Q , 

V • C + c " 2 dz/dt = tp 

in Q, 

z(; 0) = 0, C(-,0) = 0 

on f 2 , 

z = 0 

on r D x 7, 

C ' n n = 0 

on Tjv x I, 

-z - C • nn = 0 

on V R x 7, 


with data ip G L 2 (Q), $ G L 2 (Q) n . In the next proposition 
we assume that there exists a positive constant C/tab • depend¬ 
ing on the domain Q, the mesh 7 / (thus on 7 ) and the parame¬ 
ters c, t?, a, fd, S, such that for all ip G L 2 (Q), G L 2 (Q) n the 
solution (z,C) of (fTOt satisfies the stability bound: 


(ID 


^ (l+7 2 )n^ U^ 


+ 2 


+ 


^ (i + 7>y 2 c 


2 


»-V2 


1-7 

C-n? 


L 2 { ^ uR D) 


LHJT x J^) n 
/3~ 1/2 : 



( C \V2 

2 

( t?\V 2 

+ 

((1 — <5)$) C n “ 

+ 

L 2 (^«) 

(fc) * 




L 2 (^«) 


< Qtab \ 


1*1 


l 2 (QY 


+ \\cil 


II L 2 (Q) 


if go = gN = 0 (or the corresponding parts dFff of the 
boundary are empty) then the 111 (w, r) 11 |dg+ norm at the right- 
hand side can be substituted by 111 (w, t) | | |dg- 

Combining coercivity and continuity, Cea lemma gives well- 
posedness and quasi-optimality of the Trefftz-DG formulation 
irrespectively of the discrete Trefftz space V(77,) chosen. 

THEOREM 3 The variational problem <[6]i admits a unique so¬ 
lution ( Vhp,o- h p ) G V(Tft). It satisfies: 


Conditions under which bound (TTTI) holds are given in 35.51 

PROPOSITION 4 If bound ( I I II ) is satisfied, then for all Trefftz 
fields (w,r) GTffi) 


( 12 ) 


( || -1 112 

I 11 c w 


L 2 (Q) 


II L 2 (Qr 


1/2 


< C' s tab|||(w,'r)||| D G- 


Proof. We first prove the vanishing of certain jumps across 
mesh faces for the solution (z, £) of the inhomogeneous aux¬ 
iliary problem (ITOl) : \z\ t and [</] t on ,F/ pace and of [z]n and 
[C]n on J)]™. Given a hyperplane II = {nfj-x+nn f = Cn}, 
denote the scalar jump of a piecewise continuous function / as 
[/]n := /|{n*-x+n‘ I *>C , n} ^/|{n“ x+n* 1 t<Cn}- From £Kj]>, 
the fields (£, c~ 2 z) and (zej,£j), 1 < j < n (ej denoting the 
standard basis elements of M”), are in U(div Xi 4 , Q), thus their 
normal jumps vanish across any space-time Lipschitz interface 
in Q and in particular across II: 


IC' n n+ c zn nJn — [M n n)i+Cj n n]n — 0 


1 < / < n. 


In the next sections we control the L 2 (Q) norm of the error 
with a duality-type argument inspired by j 8 j Theorem 3.1], 
based on an auxiliary problem and on energy estimates. 

5.3 Auxiliary problem and error bounds in L\Q ) 

Consider the auxiliary inhomogeneous IB VP 


Thus, on time-like faces n f u = 0, the jump of z and the normal 
jump of £ vanish. On constant-time faces (n^ = 0, = ±1) 

all jumps vanish. On other planes, nfj f 0 and |njj| ^ 0, thus 

(-c 2 /<)[C • n?j]n = Mn = (-</|n£| 2 )[< • n£] n . 


I f </ 


|n^| 7 ^ c then we have immediately W n = [C-n£] n = 
0 and from above, [Cj]n = 0 for all 1 < j < n. Assumption 
(0 guarantees that njj/lnfjl > c on J 7 / 11306 , so we conclude that 
all jumps vanish. (We have simply shown that the discontinu¬ 
ities of solutions of the first-order wave equation with source 
term in L 2 (Q) n+1 propagate along characteristics.) 

Since the L 2 {Q) norm of ( w , r) G T(7/) can be computed as 


(13) ( 


-1 


w\ 


L 2 (Q) 


1 2 U /2 

Il 2 (Q) ) 
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sup 


ffo (wip + r • <f») cLcdt 


(*^)€L a W)»+ 1 (||^||22 (Q) n + \\C1P\\ L 2 {Q) ) 


1/2 1 


we now take the scalar product of (w, r) with the source terms 
(//>, 4>) of problem (1 1 01 1 and integrate by parts in each element: 


f (wijj + t ■ C) dV 

Jq 


{To) 


= E 

KdT h 


— M w dz 

, K A v ' c + ?at 


T-V: + r I dr 

Ot 


[ (wC ■ n K + T ■ n x K z- 

dcr JdK v 


wzn' 


K 


t • £ri 


Ker h 


/ F p 


t K )dS 


[wC + T2-]n + [c W2 + T-Clt dS 

' ->v--—' 

<c- 1 IMt|(7lCI+c- 1 | z |)+IH t |(7c- 1 ld+ICI) 

[ (J^+r. C 


>n " c 


V C2 Q 


dS 


=o 


[wC + tz)n d S 


= MN-C+['r] N i; 

+ f (w C ‘ n o +T-nJj ^ ) d5 

=0 on J 7 ^ 

2 




=0 on 


<[||(«;,t)||| dg ■ ^2 ( 

'(1 + 7 2 KfV /2 


(l + 7 2 )n^^/2£ 

c 




1-7 

2 


i2 ( JF space) 


+ 2 lie 1 2 1 ' 2 


i 2 (^) 


+ 2 HCII 


L 2 (^ pa “) 


+ 

a- 1/2 C • n : ' F 

2 

+ 

"Co 

1 

to 



L 2 (^ me U^ r f) 



L 2 (^r e u^) 


\ 1/2 

) C ' n r 


. (1 - <5)i?/ 


L 2 (^ 


1?\l/2 


© 


L 2 (r?) 


) 


1/2 


GD 


<© tab |||KT)||| DG (||$||i 2(Q) „ 


IIl 2 (Q) 


1/2 


Inserting this bound in the expansion (IT3l> of the L 2 (Q) norm 
of (w, r), we obtain the assertion dT2l> . □ 

From Proposition [4] and Theorem [3] it follows that the L 2 (Q) 
norm of the Trefftz-DG error is controlled by the ||| • |||dg+ 
norm of the best-approximation error: if bound dTTb is verified. 



In 0, bound (im > was proven in one space dimension on 
meshes made of rectangular elements aligned to the space- 
time axes and C s t a b was computed explicitly. Two proofs were 
given. One of them (Appendix A of 0) relies on the use of the 
exact value of ( 2 , £) in Q computed with Duhamel’s principle, 
and can not be easily extended to general domains in higher di¬ 
mensions as it require a suitable periodic extension of the IB VP 
( ITOb to R ra x (0, T). The second proof (Lemma 4.9 of 0) uses 
an energy argument to control the traces on space-like faces 
in CD and an integration by parts trick to bound the traces 
on time-like faces. In higher dimensions, the energy argument 
carries over, while the traces on time-like faces are harder to 
control. In 35.51 we prove the stability estimate (fTTT i in any di¬ 
mension, under two additional assumptions to get around the 
need to control traces on time-like faces. We first discuss in the 
next section some energy identities which will be useful later. 

5.4 Energy identities and estimates 

We call “space-like interface” a graph hypersurface 

£ = {(x,/e(x)) : X £ 0} C Q 

where /s : f2 —>■ [0, T] is a Lipschitz-continuous function 
whose Lipschitz constant is smaller than 1/c. Each space¬ 
like mesh face in © pilcc is subset of a space-like interface E. 
The future-pointing unit normal vector is debited almost ev¬ 
erywhere on E and denoted by (nf>, n\d). 

For sufficiently smooth scalar and vector belds (w, r), dehne 
their energy on a space-like interface as 

£(E ;w,t) := j (wt • n| + + M 2 )n^ dS. 

The energy on constant-time, or “flat”, space-like interfaces are 
denoted by £ (t; w. r) := £(fl x {/}; w, r), for 0 < t < T. 

For two space-like interfaces Ei, E 2 with in Q, we 

denote the volume between them and its lateral boundary as 

OEi.Ea := {(x,f) : x £ fi, / Sl (x) < t < /e 2 (x)}, 
r S i,s 2 := {(x,f) : x £ dfl, / El (x) < t < /s 2 (x)}. 

For such interfaces, the following energy identity can be veri¬ 
fied integrating by parts: 


(15) 


wt ■ n^dS 


£(E 2 ;u>,t) = £(Ei;u>,t) - J 

+ /«*.„, (( v " + ©' T+ ( v ' T+ 1 dK 
If (v, cr) is the solution of IB VP (0, then we have 


£(E 2 ;t>,er) = £(Y, 1 ;v,cr) 



v cr ■ dS 1 . 
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If <?£> = <7 jv = 9r = 0 in their domains of definition, since 
0 > 0, we have £ (£2; v, cr) < £ (Ei; v, cr), i.e. energy is dissi¬ 
pated. If moreover = 0, then £ (£2; v, cr) = £(Ei; u, er), 
i.e. energy is preserved. 

Expanding f(v hp ,a hp ) > |||(v/ l p,cr fl p)|||| )G as in (5J eq. (28)], 
it is easy to prove that the method <(6ji is dissipative and that en¬ 
ergy is dissipated by the discrete solution jumps across mesh 
interfaces: if go = gN = gn = 0, the energy of the discrete 
solution (vhp, cr/ip ) of the Trefftz-DG formulation satisfies 

£{T\v hpi (r hp ) < £(0-v o ,a o ). 

Using |nf.| < 7 on and the weighted Cauchy- 

Schwarz inequality, we have lower and upper bounds for the 
energy: if E C J^ pace u jzo y jrT j s a S p ace -iik e interface 
composed by element faces, then for all 0,1-) e L 2 (E) 1+rl 

(16) ^ J (1 + |r| 2 ) dS < £(S;w,r) 

-^/ E (1 + 7)4 & + |T|2 ) d5 ‘ 

5.5 Stability of the auxiliary problem 

We are now ready to prove bound dTTb which guarantees the 
error bound in L 2 (Q ) norm (1 1 41 . We note that, under assump¬ 
tion (ii) below, the wave speed c must be constant throughout 
Q. A mesh satisfying this assumption is depicted in FigureQ] 

PROPOSITION 5 Consider the 1BVP (ITOl) and assume that 


(i) I'd = I’.y = 0, i.e. only Robin boundary conditions are 
allowed, dfl = Tr; and 

(ii) (FI)'" 0 = 0, i.e. no time-like mesh interfaces are allowed. 


Then the stability estimate (ITU holds true with 


C; 


2 

stab 


=2 t(j 


N 




4(1+ 7 2 ) 


(1-7) 2 


+ 


r e ) 


6(1 - 5 ) 




where N is the minimal number of space-like interfaces 
Ei,..., Eat such thatT s h pace C Ui<,<iv-i and 0 < /si < 
• • • < h N -i < /sjv = T. 


GD 

< 


< 2 T 


2 J (^£(0;zX)+J J (* • C + A) dsdx^j dt 


(r 1 


=0 

II 2 

IIl 2 (Q) 


IICII 


1/2 


l 2 (QY‘ 


1 i 2 (Q) 


1*1 


L 2 (QY 


) 

\ 1/2 
J ’ 


from which 


(18) 

II c -^IIL( Q )+IICII^(q^<4t 2 (||^|| 2 l2(q) +ii$ii 2 l2(q) „). 

For every space-like mesh interface E C Jj, we control the 
corresponding term in (fTTb with the energy term: 


(19) 


^ ( 1 + 7 2 )nh y r * z +2 ^+T)n^t y/^ 


1/2; 


< 


7 

4(1+ 7 2 ) 


(1-7) 2 




L 2 ( S) 

1 
2 


' (I+7 2 )?' 


1/2 


L 2 (S) rl 


J 7) n s(^2 + IC| 2 ) dS 


~.c~ 


m 

< Cry £(E; Z, C). 


We partition the faces in J7 pace into (N — 1) interfaces J7P ace = 

Ui<j<jv-i such that 0 < /si < • • < /s w _! < T and 

denote E^r = fix {T}. We now control all terms on the space¬ 
like faces: 




1-7 


c 

„2\„t 


^ (l + 7 2 )< j 1/2 c 


L 2 (^ plce U Jff) 

2 


L*(J™**U J^) r 


til A 

< C^£(Zy,z,C) 

3=1 

UD " f . , 

<c^J2 (*•<: +a) du 

j —1 -fQnxlol.Sj 

< 2C 7 WT( ||C^|| 2 2 (Q) + ||*||i2(Q)n )• 


Proof. Applying the energy identity (fT5l > to the solution (0, (j) 
of the IBVP ( ITOl ). we have that for any two space-like interfaces 
Ei, E 2 with /ej < /e 2 , 


(17) 


f (E 2 ; z, C) < f (E i; z, C) + [ (* • C + 4>z) dV 

J Qn 1 ,X)2 


(equality holds if I'e, ,e 2 H J-j f has vanishing n-dimensional 
measure). This implies a bound in space-time L 2 norm: 


I -1 II 2 

l c Z W LHQ) 


IICII 


L 2 (QY 


= 2 


£(t-z,C)dt 


From assumptions (i) and (ii), in (ITU there are no terms on 
jziune j:D anc j jrN, so we are now ] e f t w ith the terms on 
using the Robin boundary condition —z = C ' n G , the energy 
identity (IT5l) and the L 2 (Q) stability bound (IT8l) . we have 


( C \ !/ 2 

2 

/ i) \V2 

((1-5)1?) C n ° 

+ 

L 2 (S?) 

(&) 2 




< 


1 


5(1-5) 


A? 


^ R (5|(C-n^) 2 + (l-5)^ 2 ) 


0 


= :C,5 
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A. Moiola 


zC-n^dS 

^ C s f g(0 ;z, C )- £{T;z, C) +[ (*-C + '<l’z)dV^ 

=o >o 

— Cg2T( ||C0|| £ 2(Q) + ll^lll,2(Q)n )• 

Combining this inequality with the previous one we obtain the 
assertion with Cg tab = 2T(C 1 N + Cs). □ 

Assumption (ii) in Proposition 0 requires that all the internal 
mesh faces are space-like; Figure [T] shows a mesh of this kind. 
The meshes that satisfy this condition allow the Trefftz-DG 
method to be treated as a “semi-explicit” scheme as in BJEJ: 
if the elements are suitably ordered, the discrete solution can 
be computed sequentially solving a local problem for each el¬ 
ement. This also allows a high degree of parallelism. If the 
“ tent-pitching ” algorithm of El is used to construct the mesh 
and the “macro elements” of 0 are taken as elements, then the 
mesh obtained satisfies the assumptions of Proposition 0 The 
fact that the elements obtained in this way do not have simple 
shapes such as (n + l)-simplices is not a computational diffi¬ 
culty: all integrals in the formulation © are defined on mesh 
faces, which are n-simplices, thus no quadrature on compli¬ 
cated shapes needs to be performed. This is due to the Trefftz 
property, so this advantage is not available to discretisations 
employing standard (non-Trefftz) local spaces. 



ternal mesh faces are space-like. Not all mesh elements are 2- 
simplices (triangles), but the faces are 1-simplices (segments), 
so all integrals in © are easy to compute. For this mesh, the 
parameter N in the proof of Proposition 0 is equal to 16. 

6 Polynomial Trefftz spaces 

So far we have not specified any discrete space V(7/i ): the only 
condition we imposed is the Trefftz property, i.e. V(Th) C 
T(7/i). Given p £ N, a simple choice is to define Y(Tk) := 
]\ KeTh V P (K), where Y P {K) := T(K) n P p (R n+1 ) 1+ " is 
the space of the solutions (w, t) of the wave equation in K that 
are polynomials of degree at most p. 

For high polynomial degree p, the Trefftz space V p ( K) has 
much lower dimension than the space-time, vector-valued, full 


polynomial space P P (R" +1 ) 1+ ": 

dim V r (K) = (” + ») - 1 = 0 _<p"), 

dimPP(R" +1 ) 1+ " = ( P l ") <P + ” + 1) = Op^=o(p" +1 )- 

In dimension n = 1, it has been proved in © §5.3] that these 
two spaces have the same orders of approximation in the mesh- 
width }ik := diam K and in the polynomial degree p. This is 
true both for solutions with bounded Sobolev regularity (alge¬ 
braic orders in Hk and p) and for analytic solutions (exponen¬ 
tial orders in p). Thus the order of convergence in terms of 
numbers of degrees of freedom can be considerably higher for 
the Trefftz-DG method compared to similar DG schemes based 
on standard spaces (see e.g. 13 Fig. 2]). Similar results have 
been observed in higher dimensions mm, but the proof of 
the orders of convergence is open. 

A basis of Y p (K) can be constructed using “transport polyno¬ 
mials”, or “polynomial waves”, in the form (J \Pe,j, —VPej) 
with P(j (x, t ) = (x • dj — ctY for 0 < i < p, where the 
propagation directions e R n , \d :j \ = 1, are suitably cho¬ 

sen. Non-polynomial Trefftz basis functions can easily be con¬ 
structed, for example, from /’('x-d, ct) for any smooth func¬ 
tion /:!->!. 
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